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We consider the nonequilibrium quantum vibrations of a molecule clamped between two macro- 
scopic leads in a current-carrying state at finite voltages. Our approach is based on the nonequilib- 
rium Green function technique and the self-consistent Born approximation. Kinetic equations for 
the average populations of electrons and vibrons are formulated in the weak electron-vibron coupling 
case and self-consistent solutions are obtained. The effects of vibron emission and vibronic instabil- 
ity are demonstrated using few-orbital models. The importance of the electron-vibron resonance is 
shown. 



During the past several years, nonequilibrium quantum 
transport in nanostructures and, in particular, transport 
through single molecules, has been in the focus of both 
experimental and theoretical investigations because of 
possible electronic device applications. Recently, the in- 
teraction of electrons with molecular vibrations attracted 
attention after experiments on inelastic electron trans- 
port through single molecules |1 ll 0, II II 0JiLNew 
theoretical treatments we re p resented in Refs. [d liol llll 
111 111 13 III III 13 111 III ii EH. In this Letter, we 
consider a quantum theory of nonequilibrium vibronic 
excitation. 

Basically there are two main nonequilibrium effects: 
the electronic spectrum modification |2l| and excitation 
of vibrons (quantum vibrations). In the weak electron- 
vibron coupling case the spectrum modification is usually 
small (which is dependent, however, on the vibron dis- 
sipation rate, temperature, etc.) and the main possible 
nonequilibrium effect is the excitation of vibrons at fi- 
nite voltages. We develop an analytical theory for this 
case. This theory is based on the self-consistent Born 
approximation (SCBA), which allows to take easily into 
account and calculate nonequilibrium distribution func- 
tions of electrons and vibrons. 

If the mechanical degrees of freedom are coupled 
strongly to the environment (dissipative vibron), then 
the dissipation of molecular vibrations is determined by 
the environment. However, if the coupling of vibrations 
to the leads is weak, we should consider the case when 
the vibrations are excited by the current flowing through 
a molecule, and the dissipation of vibrations is also de- 
termined essentially by the coupling to the electrons. In 
this Letter, we show that the effects of vibron emission 
and vibronic instability are important especially in the 
case of electron-vibron resonance. 

We describe a molecule coupled to free conduction elec- 
trons in the leads by a usual tunneling Hamiltonian. Fur- 
thermore, the electrons are coupled to vibrational modes. 



We do not consider Coulomb interaction to avoid fur- 
ther effects, such as Coulomb blockade and Kondo effect, 
which could dominate over the physics which we want 
to address, however self-consistent mean-field effects can 
be included easily in our approach. The full Hamilto- 
nian is the sum of the molecular Hamiltonian Hm, the 
Hamiltonians of the leads -fffl(L), the tunneling Hamil- 
tonian Ht describing molecule-lead coupling, the vibron 
Hamiltonian Hy including electron-vibron coupling and 
coupling of vibrations to the environment 

H = Hm + Hl+Hr + Ht + Hv. (1) 

A molecule (as well as a system of small quantum dots) 
is described by a set of localized states |a) with energies 
ta (tight-binding model) by the following model Hamil- 
tonian: 

= E + 6"^" W) + t^pdUp, (2) 

where d^^jda are creation and annihilation operators in 
the states \a), and (pa{t) is the (self-consistent) electrical 
potential. 

The Hamiltonians of the right (R) and left (L) leads 
are 

Hi^L(R) = ^(.fika + eLpi{t))c\f,^C,ka, (3) 
ka 
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Figure 1: (Color online) Schematic picture of the considered 
system. 
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ipi{t) are the electrical potentials of the leads, and the 
tunneling Hamiltonian 



i—L,R ka.a. 



(4) 



describes hopping between the leads and the molecule. 
Direct hopping between two leads is neglected (weak 
molecule- lead coupling case). 

Vibrations and the electron-vibron coupling are de- 
scribed by the Hamiltonian 

= ^fi^,ata, + ^^M«^(a, + aJ)4d;3. (5) 

q ap q 

Here vibrations are considered as localized phonons and 
q is an index labelling them, not the wave-vector. We in- 
clude both diagonal coupling, which describes a change of 
the electrostatic energy with the distance between atoms, 
and the off-diagonal coupling, which describes the de- 
pendence of the matrix elements tafi over the distance 
between atoms. 

We use the nonequilibrium Green function (NGF) 
method |22t , which now is a standard approach in 
mesoscopic physics and molecular electronics Q. We 
follow the formulation of Meir, Wingreen, and Jauho 
[2^. I25I , which has been already applied to the case 
of self-consistency in Ref. |2]| . 

The current in the left (z = L) or right (i = R) contact 
to the molecule is described by the well-known expression 



J^=L.B = j J ^Tr{r,(e-e(^,)(G<(e) 



(6) 



where f^{e) is the equilibrium Fermi distribution function 
with chemical potential fii, and the level- width function 
is 

ka 

The matrix lesser (retarded, advanced) Green func- 
tions of a nonequilibrium molecule G"^'-^'"^^ = G^^^'"^' 
can be found from the Dyson-Keldysh equations in the 
integral form or from the corresponding equations in the 
differential form (and references therein). 

In the standard self-consistent Born approximation, 
using the Keldysh technique, one obtains for the vibronic 
self-energies [HlIlIllElIilllllillll 

q ^ 

+M«Gf_,M«D« - 2i?« ^oM«Tr [G<M«] ) , (7) 
S<(^)(.) = ^J2J ^M^G<_M^i?< , (8) 



where G^ = 2G< + G^ - G^ is the Keldysh Green 
function, and = M^^. 

In our model the retarded vibron function is calculated 
from the Dyson-Keldysh equation 



2uj„ 



(9) 



where Il{q,uj) is the polarization operator (boson self- 
energy). The equation for the lesser function (quantum 
kinetic equation in the integral form) is 



y'-'-qul '-'-qui /-^ qui V-'^quj ^qulJ'-'-qul 



(10) 



this equation in the stationary case considered here is 
algebraic in the frequency domain. 

The polarization operator is the sum of two parts, en- 
vironmental and electronic: H^^< = H^cj^'°'™'' +H^^'°'-'. 

The environmental equilibrium part of the polarization 
operator can be approximated by the simple expressions 

H«(<=-)(g,t.) = -^7gSign(c.), (11) 

H<("™)(q,c.) - -J7,/j^(c^)sign(c.), (12) 

where 7^ is the vibronic dissipation rate, and fg{uj) is 
the equilibrium Bose-Einstein distribution function. 

The electronic contribution to the polarization opera- 
tor within the SCBA is 



H^«=»(,,c.) 



de 
2^ 



Tr(M«G<M«G:: 



+M9Gf M«G<_„) , 



H<(°')(g,tj) = -i J ^Tr (M''G<M«G>_, 



(13) 



(14) 



We obtained the full set of equations, which can 
be used for numerical calculations. We simplify these 
equations and obtain some analytical results in the vi- 
bronic quasiparticle approximation, which assumes weak 
electron-vibron coupling limit and weak external dissipa- 
tion of vibrons: 



In =lq 



2ImH^(tjg) « ( 



(15) 



So that the spectral function of vibrons can be approxi- 
mated by the Dirac (5, and the lesser function reads 

D<{q,Uj) = -~2m [{Nq + 1)5{UJ + Uq) + Nq5{uj - Uq)] , 

... 

where Nq is (nonequilibrium) number of vibrations in the 
q-th mode. So, in this approximation the spectrum modi- 
fication of vibrons is not taken into account, but the pos- 
sible excitation of vibrations is described by the nonequi- 
librium Nq. The dissipation of vibrons is neglected in the 
spectral function, but is taken into account later in the 
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kinetic equation for Nq. A similar approach to the single- 
level problem was considered recently in 0, 0, Il9| . 
The more general case with broadened equilibrium vi- 
bron spectral function seems to be not very interesting, 
because in this case vibrons are not excited. Neverthe- 
less, in the numerical calculation it can be easy taken 
into consideration. 

From the general quantum kinetic equation for vibrons 
()10|l we obtain in this limit 



+++++++ 



Nr, 



7gA^°-Imn<K) 
7, -2Imn^K) ■ 



(17) 



This expression describes the number of vibrons iV^ 
in a nonequilibrium state, — fg{LOq) is the equi- 
librium number of vibrons. In the linear approxima- 
tion the polarization operator is independent of Nq 
and — 2Imn-^((jjq) describes additional dissipation. Note 
that in equilibrium Nq = N^ because Imn<(a;g) = 
2lmn^{ujq)f'^{ujq). See also detailed discussion of vibron 
emission and absorption rates in Refs. [l^ . 

For weak electron-vibron coupling the number of vi- 
brons is close to equilibrium and is changed because of vi- 
bron emission by nonequilibrium electrons, Nq is roughly 
proportional to the number of such electrons, and the 
distribution function of nonequilibrium electrons is not 
change essentially by the interaction with vibrons (per- 
turbation theory can be used). The situation changes, 
however, if nonequilibrium dissipation —2lmll^{ujq) is 
negative. In this case the number of vibrons can be essen- 
tially larger than in the equilibrium case {vibronic insta- 
bility)^ and the change of electron distribution function 
should be taken into account self-consistently. 

In the stationary state the nonlinear dissipation rate 



7g =7q 



2Imn^(Ljg) 



(18) 



is positive, but the nonequilibrium contribution to dissi- 
pation — 2Imn^(cL'g) remains negative. 

Additionally to the vibronic quasiparticle approxima- 
tion, the electronic quasiparticle approximation can be 
used when the coupling to the leads is weak. In this 
case the lesser function can be parameterized through 
the number of electrons in the eigenstates of the non- 
interacting molecular Hamiltonian 
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Figure 2: (Color online) Vibronic emission in the symmetric 
multilevel model: voltage-current curve, differential conduc- 
tance, and the number of excited vibrons in the off-resonant 
(triangles) and resonant (crosses) cases (details see in the 
text). 



Note that in the calculation of the self-energies and po- 
larization operators we can not use ^-approximation for 
the spectral function (this is too rough and results in the 
absence of interaction out of the exact electron-vibron 
resonance). So that in the calculation we use actually 
(|19|l with broadened equilibrium spectral function. This 
approximation can be systematically improved by includ- 
ing nonequilibrium corrections to the spectral function, 
which are important near the resonance. It is important 
to comment that for stronger electron-vibron coupling 
vibronic side-bands are observed in the spectral function 
and voltage-current curves at energies ± rnoq, we do 
not consider these effects in the rest of our paper and 
concentrate on resonance effects. 

After correspondingly calculations we obtain finally 



Nr, 



N° 
q 



J2r,S Kr,5{uJq)Fr,{Fs - 1) 



(21) 



where coefficients n^jg are determined by the spectral 
function and electron-vibron coupling in the diagonal 
representation 



K,jg{ujq) 



^M^^^Aggie-iVqWlJ^^ie), 



(22) 



(19) 



we introduce the unitary matrix S, which transfer the 

r(0) 



Hamiltonian H = into the diagonal form H = 



S ^HS, so that the spectral function of this diagonal 
Hamiltonian is 



Agr,{e) = 27r(5(e - eg)Sgr,, 
where eg are the eigenenergies. 



(20) 



F„ 







'Qg'FgNq- 


Hts'Fsil+Nq) 


^Lrjri+^B.rrn'^^qri 


C7(i 


-Fg + Nq) 


+ C^s\Fg+Nq)_ 



(:t! = M^^,Agg{l,±LUq)Ml^, 



(23)- 



(24) 



here F^^^ and ff^ are the level width matrix in the diago- 
nal representation and Fermi function at energy — e(pi. 
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Figure 3: (Color online) Vibronic instability in an asymmetric 
multilevel model: voltage-current curve, differential conduc- 
tance, and the number of excited vibrons (crosses). Dashed 
line show the voltage-current curve without vibrons (details 
see in the text). 



Figure 4: (Color online) Floating level resonance: voltage- 
current curve and the number of excited vibrons (crosses). 
Dashed line show the voltage-current curve without vibrons 
(details see in the text). 



These kinetic equations are similar to the usual golden 
rule equations, but are more general. 

Now let us consider several examples of vibron emission 
and vibronic instability. 

(i) Vibron emission. First we consider the most sim- 
ple case, when the instability is not possible and only 
vibron emission takes place. This corresponds to a neg- 
ative imaginary part of the electronic polarization oper- 
ator: Imn^)(a;,) < 0. From the Eq. ^ one can see 
that for any two levels with the energies e,, > es the coef- 
ficient Kjjs is larger than k^^, because the spectral func- 
tion ^^^(e) has a maximum at e = e^-. The contribution 
of Kr,s{oJq){F^ ~ Fs) is negative if < Fg. This takes 
place in equilibrium, and in nonequilibrium for transport 
through symmetric molecules, when higher energy levels 
arc populated after lower levels. The example of such 
a system is shown in Fig.[21 Here wc consider a sim- 
ple three-level system (ei = 1. £2 = 2, es = 3) coupled 
symmetrically to the leads (Fi,, = F;?,, = 0.01). The 
current-voltage curve is the same with and without vi- 
brations in the case of symmetrical coupling to the leads 
and in the weak clcctron-vibron coupling limit (if we ne- 
glect change of the spectral function). The figure shows 
how vibrons are excited, the number of vibrons Ny in the 
mode with frequency ujq is presented in two cases. In the 
off-resonant case (green triangles) Ny is very small com- 
paring with the resonant case {ujo = £2 — ei, red crosses, 
the vertical scale is changed for the off-resonant points). 
In fact, if the number of vibrons is very large, the spec- 
tral function and voltage-current curve are changed. We 
shall consider this in a separate publication. 

(ii) Vibronic instability. Now let us consider the situ- 
ation when the imaginary part of the electronic polariza- 
tion operator can be positive: Imn''^(a;g) > 0. Above we 



considered the normal case when the population of higher 
energy levels is smaller than lower levels. The opposite 
case F2 > Fi is known as inversion in laser physics. Such 
a state is unstable if the total dissipation 7* (|18|l is neg- 
ative, which is possible only in the nonstationary case. 
As a result of the instability, a large number of vibrons is 
excited, and in the stationary state 7* is positive. This 
effect can be observed for transport through asymmetric 
molecules, when higher energy levels are populated be- 
fore lower levels. The example of a such system is shown 
in Fig.|21 It is the same three-level system as before, 
but the first and second levels are coupled not symmet- 
rically to the leads (F^i = 0.001, Tm = 0.1, Fl2 = 0.1, 
r_R2 = 0.001). The vibron couple resonantly these lev- 
els {ujq = £2 — El)- The result is qualitatively different 
from the symmetrical case. The voltage-current curve is 
now asymmetric, a large step corresponds to the resonant 
level with inverted population. 

Note the importance of the off-diagonal clcctron-vibron 
coupling for the resonant effects. If the matrix M in the 
eigen-state representation is diagonal, there is no reso- 
nant coupling between different electronic states. 

(Hi) Floating-level resonance. Finally, let us consider 
the important case, when initially symmetric molecule 
becomes asymmetric when the external voltage is ap- 
plied. The reason for such asymmetry is simply that 
in the external electric field left and right atoms feel dif- 
ferent electrical potentials and the position of the levels 

= Co''' + sfa is changed (float) with the external volt- 
age. The example of a such system is shown in Fig.0] 
Here we consider a two-level system, one level is coupled 
electrostatically to the left lead ii oc (pL, the other level 
to the right lead £2 oc tpR, the tunneling coupling to the 
leads also is not symmetrical (F^i = 0.1, Tm = 0.001, 
Tl2 ~ 0.001, Ffl2 = 0.1). The frequency of the vibration. 
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coupling these two states, is luq = 1. When we sweep the 
voltage, a peak in the voltage-current curve is observed 
when the energy difference £1—62 oc eV is going through 
the resonance ei — e2 ~ wq- 

In conclusion, we considered the excitations of quan- 
tum molecular vibrations in the nonequilibrium state and 
their influence on the voltage-current curves of a single 
molecule placed between two equilibrium leads. The im- 
portance of vibron emission and vibronic instability in 
molecular transport is demonstrated. 

We thank J. Keller and K. Richter for valuable dis- 
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Foundation under grant 1/78 340 and by the EU under 
contract IST-2001-38951. 
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